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Abstract. In this paper, we describe a geometric setting for higher-order la- 
grangian problems on Lie groups. Using left-triviahzation of the higher-order 
tangent bundle of a Lie group and an adaptation of the classical Skinner-Rusk 
' I formalism, we deduce an intrinsic framework for this type of dynamical sys- 

tems. Interesting applications as, for instance, a geometric derivation of the 
higher-order Euler-Poincare equations, optimal control of underactuated control 
systems whose configuration space is a Lie group are shown, among others, along 
Oh the paper. 

< 

I I 1. Introduction 

In 1901 [21], H. Poincare variationally deduced the equations of motion of a 
fH mechanical system specified by a lagrangian / : g — )■ M where q is the Lie algebra 

of a Lie group G. These equations are actually known as the Euler-Poincare 
^ equations. There mainly appear as the reduction of a lagrangian L : TG — ?■ M 

I — I being left or right-invariant. This procedure is called Euler-Poincare reduction 

I (see [inilSn])- The Euler-Poincare equations are 

where ac?| is the dual operator of the adjoint endomorphism ad^ '■ & Q defined by 

^ ad^^ = [C,,^] where G g. Typically, it is used the functional derivative notation 

^ 51/5^ since the equations are also valid for infinite dimensional Lie algebras. As 

^ particular examples, these equations include the equations for rigid bodies and 

J> fluids, but in the latter case, one must use infinite dimensional Lie algebras. In 

k>( the variational deduction of the Equations ([T]), it is necessary a careful analysis 

^ of the admissible infinitesimal variations to deduce the equations. More precisely, 

^ the variations are obtained by a reduction procedure of the admissible variations 

of the unreduced lagrangian L : TG — )■ M. Of course, using the corresponding 
Legendre transformation, it is possible to rewrite Equations ([T| as the Lie-Poisson 
equations on q*. Observe that, as an essential feature. Equations ([T]) involve half 
of the degrees of freedom as compared with the usual Euler-Lagrange equations 
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for L : TG M and, moreover, Equations ([T]) are first-order differential equations 
while the standard Euler-Lagrange equations are second-order ones. 

Very recently, and from different motivations, there appears a considerable inter- 
est on the extension of Equation ([T]) to higher-order mechanics (see, for instance, 
[8],[l9] as typical references for higher-order mechanics on tangent bundles). Our 
main objective in the paper is to characterize geometrically the equations of motion 
for an optimal control problem of a possibly underactuated mechanical system. In 
this last system, the trajectories are "parameterized" by the admissible controls 
and the necessary conditions for extremals in the optimal control problem are ex- 
pressed using a "pseudo-hamiltonian formulation" based on the Pontryaguin max- 
imun principle or an appropriate variational setting using some smoothness con- 
ditions [1] . Many of the concrete examples under study have additional geometric 
properties, as for instance, that the configuration space is not only a differentiable 
manifold but it also has a compatible structure of group, that is, the configuration 
space is a Lie group. In this paper, we will take advantage of this property to 
give a closed and intrinsic form of the equations of motion of our initial optimal 
control problem. For it, we will use extensively the Skinner- Rusk formalism which 
combines simultaneously some features of the lagrangian and hamiltonian classical 
formalisms and, as we will show, it is the adequate space to study the problems 
that we want to characterize [22j . Other interesting characteristic of the problems 
under study is that, for their characterization, it is necessary to use higher-order 
mechanics, that is, the phase space on the lagrangian side has coordinates which 
specify positions, velocities and accelerations. In the applications, we will use only 
second-order tangent bundles, but since the extension to /c-order tangent bundles 
does not offer special difficulties we will develop our geometric theory in this last 
case. 

Moreover, in a recent paper [13] , the authors study from a pure variational point 
of view, invariant higher-order variational problems with the idea to analyze higher- 
order geometric fc-splines with the application to image analysis. The obtain as 
main result a version of the higher-order Euler-Poincare equations. We will show 
in our paper that our analysis is complementary to the one in [13j. Assuming that 
our initial problem is invariant we deduce geometrically the same equations than 
the authors, but they use mainly variational techniques. 

The paper is structured as follows. In Section |2| we introduce some geometric 
constructions which are used along the paper. In particular, the Euler-Arnold 
equations for a hamiltonian system defined on the cotangent bundle of a Lie group 
and their extension to higher-order cases. In Section |3]T, we define the Pontryaguin 



bundle G x kQ x kg* where we introduce the dynamics using a presymplectic 
hamiltonian formalism. We deduce the fc-order Euler-Lagrange equations and, as 
a particular example, the fc-order Euler-Poincare equations. Since the dynamics 
is presymplectic it is necessary to analyze the consistency of the dynamics using 



a constraint algorithm [HI [15]. Section 3.2 is devoted to the case of constrained 



dynamics. We show that our techniques are easily adapted to this particular 
case. As an illustration of the applicability of our setting, we analyze the case 
of underactuated control of mechanical systems in Section |4] and, as a particular 
example, a family of underactuated problems for the rigid body on S0{3). 
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2. Geometric preliminaries 

2.1. Euler- Arnold equations. (See, for instance, 0) Let G be a Lie group. 
Consider the left-multiplication on itself 

GxG — >G, {g,h) ^ £g{h) = gh . 

Obviously £g is a diffeomorphism. (The same is valid for the right-translation, but 
in the sequel we only work with the left-translation, for sake of simplicity). 

This left multiplication allows us to trivialize the tangent bundle TG and the 
cotangent bundle T*G as follows 

TG ^ GxQ, {g,g)^{g,g-'g) = {g,Tg£g-.g) = ig,0 , 
T*G ^ GxQ*, {g,ag)^ig,T:£g{ag)) = {g,a) , 

where q = TeG is the Lie algebra of G and e is the neutral element of G. In the same 
way, we have the following identifications: TTG = G x 3q, T*TG = G x g x 2q*. 
TT*G = G X s* X X 0* and T*T*G = G x 3q*. 

Using this left trivialization it is possible to write the classical hamiltonian equa- 
tions for a hamiltonian function H : T*G ^ from a different and interesting 
perspective. 

For instance, it is easy to show that the canonical structures of the cotangent 
bundle: the Liouville 1-form 6g and the canonical symplectic 2-form ojg, are now 
rewritten using this left-trivialization as follows: 

i0G)ig,a)i^l,'yi) = (2) 

with {g,a) E G x q*, where G and Ui G 0*, z = 1,2 and we have used the 
previous identifications. Observe that we are identifying the elements of TagT*G 
with the pairs (^, i^) G x 0*. 

Therefore given the hamiltonian H : T*G = G x q* — M, we compute 

j^{9,(^)] ,^2) + {'^2,-^{g,a)) , (4) 

since (c/,a) G 0** = 0. 

We now derive the Hamilton's equations which are satisfied by the integral curves 
of the Hamiltonian vector field Xh on T*G. After left-trivialization, XH{g,<y) = 
(^1, z/i) where ^1 G and Ui G 0* are elements to be determined using the Hamil- 
ton's equations 

ix„ujG = dH . 
Therefore, from expressions (|3| and ^ we deduce that 



ui = -£*g ( ^(51, a) ) + ad\^a . 
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In other words, taking g = g^i we obtain the Euler- Arnold equations: 

da da 



If the Hamiltonian is left-invariant, that, is : 0* — ?■ M where h{a) = H{e,a) then 
we deduce that: 

6h 

9 = 9^, 
da 

The last equation is known as the Lie-Poisson equations for a hamiltonian 

2.2. Higher-order tangent bundles. In this section we recall some basic facts 
of the higher-order tangent bundle theory. Along the section we will particularize 
this construction to the case when the configuration space is a Lie group G. For 
more details see [HI [IS]- 

Let Q be a differentiable manifold of dimension n. It is possible to introduce an 
equivalence relation in the set (7*^(1^, Q) of fc-differentiable curves from M to Q. By 
definition, two given curves in Q, 71 (t) and 72 (t), where t G {—a, a) with a G M 
have contact of order A; at go = 7i(0) = 72(0) if there is a local chart {ip, U) of Q 
such that qo E U and 



^(<^o7i(t))|^^^ = — (<^o72(t)) 



t=0 

k/ 



for all s = 0, k. This is a well defined equivalence relation in C (M, Q) and the 
equivalence class of a curve 7 will be denoted by [7)0'^^ The set of equivalence classes 
will be denoted by T^^^Q and it is not hard to show that it has a natural structure 

of differentiable manifold. Moreover, Tq : T^^^'Q — > Q where Tq ^[7]o'^''^ = 7(0) is 

a fiber bundle called the tangent bundle of order k of Q. 

In the case when the manifold Q has a Lie group structure, we will denote 
Q = G and we can also use the left trivialization to identify the higher-order 
tangent bundle T^^^G with G x kQ. That is, if 5? : / — > G is a curve in C^''^(]R, G): 

Jik) . rp{k)Q y GxkQ 



[g]i'^ ^ igiO),g-HO)m,i {g-\t)g{t)), ■ ■ ■ , §^ {9-\mm 



t=o 



It is clear that T'^'^^ is a diffeomorphism 

We will denote by ^{t) = g^^{t)g{t). Therefore 

r^'\[9t^) = i9,^,i...,i^'-'^), 

where 

^^'\t) = ^^ig-\tm), 0<l<k-l 

and ^(0) = ^,^('^0) = < I < k - 1. We will indistinctly use the notation 
C^"'* = ^5 ^^^^ = ^5 where there is not danger of confusion. 
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We may also define the surjective mappings r^'^^ : T^^^G — >■ T^^^G^ for Z < /c, 
given by Tq^^ ^blo*^^^ — [g\t^ ■ Witfi tfie previous identifications we have that 

4'\g{Q)mm. . . . , = (^(0), e(o), e(o), . . . , e^^-^Ho)) 

It is easy to see that T^^^G = G x g, T^^^G = G and t^^'^ = r^. 

Now, we consider the canonical immersion '■ T^^^G — )■ TiT^^'^^^G) defined as 
jk{[gfo^) = [9^''~^'^]o \ where g^'"'^'^ is the lift of the curve g to T^'^'^'^G; that is, the 
curve ^('^-i) : M -> T^^'^^^G is given by g^''-^\t) = [gt]o~^^ where gt{s) = g{t + s). 
Using the identification given by T*^^) we have that: 

Z'^) : GxkQ ^ Gx{2k-1)Q 

{g,tL...,e'-'^) ^ {g,C,L---,C^'-'^;C,i---,e'^) 

where we identify T{T^^-^^G) = T{G x {k - 1)q) = G x {2k - 1)0, in the natural 
way. 

2.3. Higher-order Euler-Arnold equations on 

T*{T^k-i)Qy Combining the 

results of the two previous subsections we have that 

= T*{G x{k- 1)0) = T*G x{k- 1)T*Q = G x {k - 1)q x kg* . 

For developing our geometric formalism for higher-order variational problems 
on Lie groups we need to equip the previous space with a symplectic structure. 
Thus, we construct a Liouville 1-form ^Gx(fc-i)0 and a canonical symplectic 2-form 
ujGx{k-i)s after the left-trivialization that we are using. Denote by £ e (A; — 1)0 
and a G kg* with components ^ = {^^^\ . . . , C*-''"^-') and a = (ckq, . . . , cxk-i)- Then, 
after a straightforward computation we deduce that 

(6'Gx(fc-l)fl)(g,4,a)(6,»^^) = 

fc-1 



1=0 



where G kg and iv" e A;0*, a = 1,2 with components = (^a )o<j<fc-i and 
= ('^("))o<i<fc-i where each component ^a'' € and u^-^ e 0*. Observe that ao 
comes from the identification T*G = G x g*. 

Given the hamiltonian H : T*T^^-^^G = Gx{k-l)gx kg* — > R, we compute 

dif(,,,,«)(6,^^) = (i^;(^(^7,^,«)),ef) + E(^(^'^'")'^^'^) 



+{u^^-^{g,^,a)) 

As in the first subsection, we can derive the Hamilton's equations which are 
satisfied by the integral curves of the Hamiltonian vector field Xh defined by 
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XH{g,^,a.) = ($1,1^^). Therefore, we deduce that 

OOL 

5H 

= --^i9,^,Oi), 0<i<k-2. 

In other words, taking g = g^^^^ we obtain the higher-order Euler- Arnold 
equations: 

9 = 9^{9,^,ol) , 

Otto 

-{g,^,a), l<t<k-l, 



dt 6ai 
dao 



ig,tct), 0<t<k~2. 



dt 

dt ^ 

3. On the GEOMETRY OF HIGHER-ORDER VARIATIONAL PROBLEMS ON LlE 

GROUPS 

In this section, we describe the main results of the paper. First, we intrinsically 
derive the equations of motion for Lagrangian systems defined on higher-order 
tangent bundles of a Lie group and finally, we will extend the results to the cases 
of variationally constrained problems. 

3.1. Unconstrained problem. In 1983, it was shown by R. Skinner and R. Rusk 
[22] that the dynamics of an autonomous classical mechanical system with la- 
grangian L : TQ M can be represented by a hamiltonian presymplectic system 
on the Whitney sum TQ®T*Q (also called the Pontryaguin bundle). The Skinner- 
Rusk formulation can be briefly summarized as follows. Denoting the projections 
of TQ © T*Q TQ and TQ © T*Q T*Q by pri and pra, respectively, we define 
the presymplectic 2- form Q = pr2UjQ (where ojq is the canonical symplectic 2-form) 
and the hamiltonian function H : TQ © T*(5 — i- M by 

H{vg, ag) = {ag, Vq) - L{vq), aq G T*Q and Vq G TqQ . 

Then, one can consider the following equation 

ix^ = dH , 

which describes completely the equations of motion of the Lagrangian system. 
Moreover, if the Lagrangian is regular then there exists a unique solution Z which is 
tangent to the graph of the Legendre map. This formalism is quite interesting since 
it allows us to analyze the case of singular lagrangians in the same framework and, 
more important, it also provides an appropriate setting for a geometric approach 
to constrained variational optimization problems. 



HIGHER-ORDER VARIATIONAL PROBLEMS ON LIE GROUPS 



3.1.1. The equations of motion. Now, wc will give an adaptation of the Skinner- 
Rusk algorithm to the case of higher-order theories on Lie groups. We use the 
identifications 

T^'^^G = GxkQ, 

j.*j.{k-l)Q ^ Gx{k-1)QX kQ* . 

Consider the higher- order Pontryaguin bundle 

Wq = T^^^G x^ik-DG T*T^^-^^G = GxkQxkQ* , 
with induced projections 

where, as usual, £ = (^^°\ . . . , ^^''~^^) e {k — 1)q and a — (ckq, . . . , Oik-i) G kQ*. 

Wo ^ G X kg x kg* 

pn ^^^^ pr2 



G X kg 



(fc-i,fc) 



Gx{k-1)QX kQ* 



Gx{k- 1)0 



For developing the Skinner and Rusk formalism it is only necessary to construct 
the presymplectic 2-form fl^^ by flwo = P''^2'^Gx{k-i)3 and the hamiltonian function 

H(g,t^^'-'\cx) = - L(g,te-'^) . 

i=0 

Therefore 

i=0 

where ^ kg, i/" G kg*, and ^i''^ G 0, a = 1,2. Observe that ^i''^ docs not 
appear on the right-hand side of the previous expression, as a consequence of the 
presymplectic character of ^Wo- Moreover, 



dH, 



'"^ SL 
5p 
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Therefore, the intrinsic equations of motion of a higher-order problem on Lie 
groups are now 

ix^Wo = dH . (5) 

If we look for a solution X{g,$,,^^'^~^\cx.) = (^i,^^'' of Equation (5) we 

deduce: 



6L 

and the constraint functions 



Observe that the coefficients C,i are still undetermined. 
An integral curve of X, that is a curve of the type 

t {g{t),at), • • • , ^^'-'Ht), ao{t), , 
must satisfy the following system of differential-algebraic equations (DAEs): 

9 = 9^, (6) 



dt 
~dt 



1<,<A;-1, (7) 

+adla,, (8) 

^ = j^)i9,te'-'')-a., 0<.<fc-2, (9) 

= jl^)^9.i^^''-'') ■ (10) 



If A; > 2, combining Equation (jlOj) with the (|9j) for i = A; — 2, we obtain 

d 5L _ 5L _ 

Proceeding successively, now with i = k — 3 and ending with z = we obtain the 
following relation: 

This last expression is also valid for k > 1. Substituting in the Equation ^ we 
finally deduce the k-order trivialized Euler- Lagrange equations: 

Of course if the Lagrangian L : T^'^^G = G x kg — > M is left-invariant, that is 
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for all g,h & G, then defining the reduced lagrangian I : kg — > M by 
we write Equations ([TI| as 

i=0 ^ 

which are the k-order Euler-Poincare equations (see, for instance, [13]). 

3.1.2. The constraint algorithm. Since f2vi/o presymplectic then ([s]) has not solu- 
tion along Wq then it is necessary to identify the unique maximal submanifold Wj 
along which ([s]) possesses tangent solutions on Wf. This final constraint subman- 
ifold Wf is detected using the Gotay-Nester- Hinds algorithm [IS]. This algorithm 
prescribes that Wf is the limit of a string of sequentially constructed constraint 
submanifolds 

... ^Wk^ ■■■ ^W2^Wi^Wo . 

where 

Wi = ^^xeGxkgxkg* \ \ i/^) = 

with i > 1 and where 

(T.W._i)^ = {(6,d'\^') e ik + l)g X kg* I ^^v^„(x)((6,d'\^'), (6, =0 

V(6,ef,^')eT.l^.-i}. 

where we are using the previously defined identifications. If this constraint algo- 
rithm stabilizes, i.e., there exists a positive integer k E N such that Wk+i = Wk 
and dimVTfc > 1, then we will have at least a well defined solution X on Wf = Wk 
such that 

From these definitions, we deduce that the first constraint submanifold Wi is 
defined by the vanishing of the constraint functions 

SL 

Applying the constraint algorithm we deduce that the following condition, if 
k>2: 

In the particular case A; = 1, we deduce the equation 

In both cases, these equations impose restrictions over the remainder coefficients 
d*^^ of the vector field X. 
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If the bilinear form c^(fc-n £(fc-i) : x g — )• M defined by 



d 



=0 ds 



s=0 



is nondegenerate, we have a special case when the constraint algorithm finishes 
at the first step Wi. More precisely, if we denote by flwi the restriction of the 
presymplectic 2-form Q to Wi, then we have the following result: 



Theorem 3.1. {Wi,Qwi) ^■^ ^ symplectic manifold if and only if 

5^L 



(13) 



is nondegenerate. 

3.2. Constrained problem. 

3.2.1. The equations of motion. The geometrical interpretation of constrained prob- 
lems determined by a submanifold M of G x fcg, with inclusion i^ : M, G x kg 
and a lagrangian function defined on it, Lj^ : M — )■ M, is an extension of the 
previous framework. First, it is necessary to note that for constrained system, in 
this paper, we understand a variational problem subject to constraints (vakonomic 
mechanics), being this analysis completely different in the case of nonholonomic 
constraints (see fTTl 17]). 

Given the pair (M, Lm) we can define the space 

Wo = Mxkg* . 

Take the inclusion i^^^ : Wq ^ G x kg x kg*, then we can construct the following 
presymplectic form 

^Wo = (^^^2 o %o)*^Gx(fc-l)0xfcfl* , 

and the function H : Wq — )■ M defined by 

k-l 
i=0 

where {g,^,^^''~^^) G M. 

With these two elements it is possible to write the following presymplectic sys- 
tem: 

txn^, = dH . (14) 
This then justifies the use of the following terminology. 

Definition 3.2. The presymplectic Hamiltonian system {Wq, ^y^^, H) will be called 
the variationally constrained Hamiltonian system. 

To characterize the equations we will adopt an "extrinsic point of view" , that is, 
we will work on the full space Wq instead of in the restricted space Wq (see next 
section for an alternative approach). Consider an arbitrary extension L : Gxkg ^ 
M of L]vt : M — 7- M. The main idea is to take into account that Equation (14) is 
equivalent to 

ix^Wo ~ dH e ann TWq , 

X e tWo , 
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where ann denotes the annihilator of a distribution and H is the function defined 
on Section ISTTl 

Assuming that M is determined by the vanishing of m-independent constraints 

then, locally,ann TWq = span {d^^} , and therefore the previous equations is 
rewritten as 

ix^Wo -dH = X^d^^ , 
X(<l>^) = 0, 

where are Lagrange multiphers to be determined. 

If X{g, ^, ck) = (^1, ^[''\ u^) then, as in the previous subsection, we obtain 

the following prescription about these coefficients: 

if = o<z<k-i, 

(5L 5^^\ 



1 5L , 5$^ 



k-2 



and the algebraic equations: 



SL (5$^ 

$^ = . 

The integral curves of X satisfy the system of differential-algebraic equations 
with additional unknowns (A^): 



9 = 9^ 



e«, i<t<k-i, 



dt 

dai+i 5L 5$^ 



^ <f ).H-g|;,.«..^d- 



SL ^ 5$^ 



$^ = . 
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As a consequence we finally obtain the k-order trivialized constrained Euler- 
Lagrange equations, 



1=0 



dp 



SL 



A 



(5$ 



Al 



5L ^ 5$^ 



(15) 



If the Lagrangian L : T^^^G = Gxkg — > M and the constraints : G x kg 
M., 1 < A < m are left-invariant then defining the reduced lagrangian I : kg — 
and the reduced constraints (p^ : kg ^ M. we write Equations (15) as 



i=0 



'dt 



dp 



61 



lAI 



A, 







(16) 



3.2.2. The constraint algorithm. As in the previous subsection it is possible to 
apply the Gotay-Nester algorithm to obtain a final constraint submanifold where 
we have at least a solution which is dynamically compatible. The algorithm is 
exactly the same but applied to the equation (14). 

Observe that the first constraint submanifold Wi is determined by the conditions 



6L 



. 



A. 



5$^ 



If we denote by flyy-^ the pullback of the presymplectic 2-form fi^y^ to Wi, it is 
easy to prove the following 



Theorem 3.3. {Wi,Q^^J is a symplectic manifold if and only if 







(17) 



is nondegenerate, considered as a bilinear form on the vector space g x 



4. An APPLICATION: UNDERACTUATED CONTROL SYSTEMS ON LlE GROUPS 

A Lagrangian control system is underactuated if the number of the control inputs 
is less than the dimension of the configuration space (see [6] and references therein). 
We assume that the controlled equations are trivialized where L : G x g ^M. 

d (5L\ ..(5L\ dL 

where we are assuming that {e"} are independent elements on g* and {ua) are the 
admissible controls. Complete it to a basis {e°',e^} of the vector space g*. Take 
its dual basis {cj} = {ca, e^} on g with bracket relations: 

The basis {cj} = {6^,6^} induces coordinates (t/",t/^) = (y*) on g, that is, if e G 
then e = y^Ci = y'^Ca + y^CA- In g*, we have induces coordinates {pa,Pa) for the 
previous fixed basis {e*} 
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In these coordinates, the equations of motion are rewritten as 
d f dL\ idL . 6L . 



dt \dy^ J '"^ dyi ' Hg 

With these equations we can study the optimal control problem that consists 
on finding trajectories {g{t),u°'{t)) of state variables and control inputs satisfying 
the previous equations from given initial and final conditions {g{to),y''{to)) and 
{g{tf),y''{tf)), respectively, and extremizing the functional 

C{g{t),y\t),u''{t))dt 

to 

Obviously (see |1],[I2] and reference therein) the proposed optimal control prob- 
lem is equivalent to a variational problem with second order constraints, deter- 
mined by the lagrangian L : G x 20 — M given, in the selected coordinates, by 

subjected to the second-order constraints 

- 1 (^) - • 

which determine the submanifold M of G x 20. 

Observe that from the constraint equations we have that 



dyAQyB^ dyAQyb^ dyj ^ 9^g^ 

Therefore, assuming that the matrix {Wab) = (a^^^) regular we can write 



the constraint equations as 



= G^{g,y\r) 
where W^^ = {Wba)'^- 



Wo = Mx 20* 




G X 



This means that we can identify TM = G x span {(cj, 0, 0), (0, Cj, 0), (0, 0, Cq)} 
where (e„ 0, 0), (0, e„ 0), (0, 0, e,) G 30. 
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Therefore, we can choose coordinates ((7,?/*,^") on M. This choice allows us to 
consider an "intrinsic point view", that is, to work directly on Wq = M x 2q* 
avoiding the use of Lagrange multipliers. 

Define the restricted lagrangian L]vt by L]vt = L \]^: M, ^ M and take induced 
coordinates on Wq are 7 = ig,y\y"',Pi,Pi)- Consider the presymplectic 2-form on 

Wo, ^Wo = (P^2 0%J*(i^Gxg). 

Using the notation (ei)o = (cj, 0, 0; 0, 0) G 30 x 2^* and, in the same way (ej)i = 
(0, a, 0; 0, 0), {ea)2 = (0, 0, e„; 0, 0), {e% = (0, 0, 0; e\ 0) and {e% = (0, 0, 0; 0, e') 
then the unique nonvanishing elements on the expression of 

(^Wo)7((ei)o, (ej)o) = GijPk, 

(Q^J,((e,)o, (6^)3) = -{n^M{e%, (e,)o) = Sf , 

Taking the dual basis (e*)o = (e*, 0, 0; 0, 0) G 30* x 20 and, in the same way (e*)i = 
(0, e\ 0; 0, 0), (e")2 = (0, e^ 0; 0, 0), (e,)3 = (0, 0, 0; a, 0) and (e,)4 = (0, 0, 0; 0, a) 
we deduce that 

(O^J = {e% A (e,)3 + (eOi A (e,)4 + ^ej,.pfe(e0o A (e^)o 

Moreover 
and, in consequence. 



The conditions for the integral curves t {g{t),y'^{t),y"'{t),pA{t),pA{t)) of a vector 
field X satisfying equations ^x^pFo ~ 

dg 



dt 

dy'^ 

~dt 
dy^ 
dt 



g{y\t)ei) (18) 

y^ (19) 
G'^(^,Z/^r) (20) 



(^n^-^^^ + (21) 



Pi + -pB^:r— (22) 



dt dy^ dy 

^Zm _~ _ d. 

Qya Qya Qya ' ^" dy^dy"" 



Pa = i^r7;r-Ps^ = ^ + ^ pbt^::^^ (23) 

(24) 
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As a consequence we obtain the following set of differential equations: 
dg 



dt 
dt 



G^{g,y\r) 




which determine completely the dynamics. 
If the matrix 



Qya.Qyb 

is regular then we can write the previous equations as a explicit system of third- 
order differential equations. It is easy to show that this regularity assumption is 
equivalent to the condition that the constrain algorithm stops at the first constraint 
submanifold Wi (see t2j,[9j, (TU] and reference therein for more details). 

4.1. Example: optimal control of an underactuated rigid body. We con- 
sider the motion of a rigid body where the configuration space is the Lie group 
G = S0{3) (see ^^). Therefore, TS0{3) ~ 50(3) xso(3), where so(3) = ig 
the Lie algebra of the Lie group SO {3). The Lagrangian function for this system 
is given by L : SO (3) x so (3) R, 

L{R, ^1,^2, ^^3) = I [ml + ml + mVj . 

Now, denote by t — )■ R{t) G 5*0(3) a curve. The columns of the matrix R{t) 
represent the directions of the principal axis of the body at time t with respect to 
some reference system. Now, we consider the following control problem. First, we 
have the reconstruction equations: 

-Vl:,{t) Vtiit) 



R(t) = R{t) 











R{t) {ni{t)Ei + n2{t)E2 + ns{t)Es 
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Eo:= \ 0, 




and the equations for the angular velocities Qi with i = 1,2, 3: 

m2{t) = {h-h)nsit)niit) + u2it) 
msit) = (/l-/2)^]l(t)^]2(^) 

where Ii,l2,h are the moments of inertia and ui,U2 denotes the applied torques 
playing the role of controls of the system. 

The optimal control problem for the rigid body consists on finding the trajecto- 
ries (-R(t), ^{t), u{t)) with fixed initial and final conditions {Rito), fi(to)), (Ritf), ^(tf)) 
respectively and minimizing the cost functional 

A = [ e{n,Ui,U2)dt = [ [ci{ul + ul) + C2{fll + nl + fll)] dt, 
Jo Jo 

with Ci,C2 > 0. The constants Ci and C2 represent weights on the cost functional. 
For instance, ci is the weight in the cost functional measuring the fuel expended 
by an attitude manoeuver of a spacecraft modeled by the rigid body and C2 is the 
weight given to penalize high angular velocities. 

This optimal control problem is equivalent to solve the following variational 
problem with constraints (|1],[I2]), 

mina= / L{n,tl)dt 
Jo 

subject to constraints I^^z ~ {h — -^2)^1^2 = 0, where 

L{n, n) = e (n, - {h - 13)^2^^, m2 - {h - /i)^^3^^i) ■ 

Thus, the submanifold M of G x 25o(3), is given by 

M = {(i?, n, n)\n3= (^^^y^^ ^1^2}. 

We consider the submanifold Wq = M x 2so*(3) with induced coordinates 

fil, 1^2,^3, ^1,^2, Pi, P2,P3, Pi, P2,P3)- 

Now, we consider the restriction given by 



— Ci 



(jifii - {I2 - 13)^2^3)' + [m2 - ih - h)^3^, 



+C2 {nl + + nl) 



For simplicity we denote by = ^^j^QiQ2- 

Then, we can write the equations of motion of the optimal control for this 
underactuated system. For simplicity, we consider the particular case Ii = I2 = 
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/s = 1 then the equations of motion of the optimal control system are: 



dt \ dt ' ' df^ dt^ 

d^pz dQ3 ^ . . d^VLi ^ . .d^Q 



2c2— ^ - 2cifi2(t)— ^ + 2cifii(t)- 



dt^ dt ' ' dt^ ' ' dt^ 

dQs 







dt 

If we consider the rigid body as a model of a spacecraft then we observe that this 
particular cost function is taking into account both, the fuel expenditure (ci) and 
also is trying to minimize the overall angular velocity (02). In Figures g and 
([2]) we compare their behavior in two particular cases: Ci = 1/2, C2 = 1/2; and 
ci = 1/2, C2 = 0. 

In all cases we additionally have the reconstruction equation 

R{t) = R{t) ini{t)Ei + n2{t)E2 + ^3(^)^3) • 
with boundary conditions R{to) and R{tf). 



The case ci = and C2 = 1, that is, we only try to minimize the overall angular 
velocity (see [23] for the fully-actuated case) is singular. We obtain the following 
system 

d^_ dih ^ „ 
dt^ dt 
dQs 

— ^ = 0. 

dt 

Observe that in this case it is not possible to impose arbitrary boundary conditions 
{R{tQ),Q{to)) and {R{tf),fl{tf)) although it is always possible to find a trajectory 
verifying initial and final attitude conditions R{to) and R{tf). 

5. Conclusions and future work 

We have defined following an intrinsic point of view the equations of motion for 
variational higher-order lagrangian problems with constraints. As a particular case, 
we obtain the higher-order Euler-Poincare equations (see As an interesting 

application we deduce the equations of motion for optimal control of underactuated 
mechanical systems defined on Lie groups. These systems appear in numerous 
engineering and scientific fields, as for instance in astrodynamics. In this sense we 
study the attitude control of a satellite modeled as a classic rigid body. 
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Figure 1. Angular velocity values for initial conditions satisfying 
fii(O) = Qi{4) = 0, i = 1,2 and fixed values of R{0) and i?(4). 

These techniques admits an easy generalization for the case of discrete systems. 
As an illustration , consider the second order tangent bundle of a Lie group left 
trivialized as T^^^G G x 20. We choose its natural discretization as three copies 
of the group (we recall that the prescribed discretization of a Lie algebra q is its 
associated Lie group G). Consequently, we develop the discrete Euler-Poincare 
equations for the discrete Lagrangians defined on : G x G x G — M. Define 
Wk = g'^^gk+i- Taking variations for Wk, where we denote = g^^^gk-, we obtain 

5Wk = -gk^Sgkgk^gk+i + gk^6gk+i 
= -T,kWk + gk^gk+igkli^gk+i 
= -EkWk + Wk^k+i, 
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Figure 2. Comparison of the functions 1/2 (fif (t) + ^ll{t) + ^l{t)) 
(left) and 1/2 {ul{t) + ul{t)) (right) in both cases 



where gk, Wk G G and S^. G 0. 

The equations of motion are the critical paths of the discrete action 

N-2 

min ^ Ld{gk, Wk, Wk+i) 

k=0 

with boundary conditions Sq = Si = Sjv-i = S^r = since we are assuming that 
9o, 9ii 9n-i and gjy fixed. Therefore, after some computations we can obtain the 
equations, 

ll^_DiLd{gk-i, Wk-i, Wk) + rw^_D2Ld{gk-i, Wk-i, Wk) 
-r\y^D2Ld{gk, Wk, Wk+i) - r*w^D^Ld{gk-i, Wk-i, Wk) 
+r^^_D^Ld{gk-2, Wk-2, Wk-i) = 

which are the discrete second-order Euler-Lagrange equations. 

Moreover, in a future paper we will generalize the presented construction of 
higher-order Euler-Lagrange equations to the case of Lie algebroids. This abstract 
approach will allows us to intrinsically derive the equations of motion for dif- 
ferent cases as, for instance, higher Euler-Poincare equations, Lagrange-Poincare 
equations and the reduction by morphisms in a unified way. We will generalize the 
notion of higher-order tangent bundle to the case of higher-order Lie algebroids (or, 
more generally, anchored bundles) using equivalence classes of admissible curves 
and extending the ideas introduced in [T7] . 

We will analyze in a future paper these and other related aspects. 
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